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ABSTRACT 

t"**- ' We have conducted a deep and uniform 1.1mm survey of the GOODS-N field with 

AzTEC on the James Clerk Maxwell Telescope (JCMT). Here we present the first 
\Q ■ results from this survey including maps, the source catalogue, and 1.1mm number- 

counts. The results presented here were obtained from a 245 arcmin 2 region with near 
uniform coverage to a depth of 0. 96-1. 16mJy beam -1 . Our robust catalogue contains 
■ 28 source candidates detected with S/N ^ 3.75, only ~l-2 of which are expected to 

be spurious detections. Of these source candidates, 8 are also detected by SCUBA at 
850 /im in regions where there is good overlap between the two surveys. The major 
advantage of our survey over that with SCUBA is the uniformity of coverage. We 
^ | calculate number counts using two different techniques: the first using a frequentist 

parameter estimation, and the second using a Bayesian method. The two sets of results 
are in good agreement. We find that the 1.1mm differential number counts are well 
described in the 2-6 mJy range by the functional form dN/dS = N' (S' / S)exp(— S/S') 
with fitted parameters S' = 1.25 ± 0.38 mJy and dN/dS = 300 ± 90mJy _1 deg -2 at 
3 mJy. 

Key words: instrumentatiomdetectors, sub-millimetre, galaxies:starburst, galax- 
ies:high redshift 
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1 INTRODUCTION 

Identifying and studying the galaxies at high redshift that 
will evolve into today's normal and massive galaxies re- 
mains a major goal of observational astrophysics. Galaxies 
discovere d in deep sub-millimetre and mm-wavelength sur- 
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are generally thought to be dominated by dusty, possibly 



merger-induced starburst systems and active galactic nuclei 
(AGN) at redshifts z > 2 with star format ion rates as high 
as SFR ~ 1000 MQyr -1 (| Blain et al.ll2002l U The high areal 
number density of these sub-mm and mm-detected galax- 
ies (SMGs), combined with their implied high star forma- 
tion rates and mea s ured FIR luminositi es (Lfir. ~ 10 12 Lq, 
iKovacs et all 120061: | Coppin et al.| [2008). makes their esti- 
mated contribution to both the global star formation den- 
sity and the su b-mm background radiation as high as 50% 
at z ~ 2 (e.g.. iBorvs et al.ll2003l ; IWall et alj|2008h . Their 
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observed number cou nts imply strong evolution between 
z = 2 and today ( e.g. IScott et al.ll2002l ; iGreve et all 12004 
ICoppin et al . 2006). The high star formation rates at early 
epochs of SMGs generally match the expectation for rapidly 
forming elliptical galaxies, a view supported by the high rate 
of mergers seen loc ally in samples of u ltra-luminous infrared 
galaxies (ULIRGs; [b ornc et alK 2000). which are plausible 
local counterparts of distant SMGs. Together, these charac- 
teristics have led many observers to surmise that SMGs are 
likely to evolve into the massive galaxies observed locally 
(e.g iDunlop et all Il994l ; ISmail et alj Il997l ; iBertoldi et all 
120071 ) and may hold important clues to the processes of 
galaxy and structure formation in general at high redshift. 

GOODS-N is one of the most intensively studied extra- 
galactic fields, with deep multi-wavelength photometric cov- 
erage from numerous ground-based and s pace-based facili- 
ties. These include Chandra in the X-rav ([Alexander et al.l 
120031 ). HST in the optical and NIR (|Giavalisco et al.ll2004 ). 
Spitzer in the NIR-MIR (Chary et al. in prep., Di ckinson et 
al. in prep.), and the Very Large Array in the radio (| Richards! 
2000, Morrison et al. in prep.), as well as highly complete 
spectroscopic surveys from ground-b ased observatories (e.g. 
IWirth et al]|2004 ICowie et al.ll2004 ). This field is therefore 
ideally suited for deep mm-wavelength studies of SMGs: the 
extensive coverage in GOODS-N allows the identification 
of SMG counterparts in X-ray, UV, optical, IR, and radio 
bands, as well as constraints on photometric redshifts and 
investigation of SMG power sources and evolution. 

Deep mm surveys of blank fields are needed in or- 
der to constrain the faint end of the SMG number counts, 
while large areal coverage is required to constrain the bright 
end. Together they provide strong constraints on evolu- 
tionary scenarios. Previous sub-mm surveys of GOODS- 
I S! have been carri e d out with SCUBA on the JCMT 
llHughes et"al1 Il998l: iBarger et al.| l200d : iBorvs et al] 120031 : 
IWang et all2004l : IPope et al.ll2005h . The 'Super-map' of the 
GOODS-N field, which was assembled from all available 
JCMT shifts covering the field, contains 40 robust sources 
at 850 /im down to an average sensitivity of 3.4 mjy (la ) 
and covers 200 arcmin 2 (IBorvs et al.ll2003l ; |Pope et alj2005t) . 
However, the r.m.s. is highly non-uniform ranging from 0.4 
mjy to 6 mjy (see Fig. [TJ . That non- uniformity presents se- 
rious complications for comparisons with multi-wavelength 
data. 

In this paper we present a new 1.1 mm survey of the 
GOODS-N field made with AzTEC l| Wilson et al .11 20081 ) at 
the 15-m James Clerk Maxwell Telescope (JCMT) on Mauna 
Kea, Hawaii. This map is the deepest blank-field survey 
undertaken during the AzTEC/ JCMT observing campaign, 
and is one of the largest, deepest, and most uniform mm- 
wavelength maps of any region of the sky. Our map cov- 
ers 245 arcmin 2 and completely encompasses the 16.5' x 10' 
Spitzer GOODS-N field and all of the previous GOODS- 
N sub-mm a nd mm-wavelength f ields, including the original 
HDF map of lHughes et~afl (| 19981 ) and the SCUBA GOODS- 
N 'Super-map' (indicated in Fig. [ T] here and presented in 
IBorvs et al]|2003l ; IPope et al.ll2005l ). The large number and 
high stability of the AzTEC bolometers has enabled us to 
produce a map with small variations in r.m.s., from 0.96- 
1.16 mjy, across the 245 min 2 field. This uniformity is a 
drastic improvement over the SCUBA GOODS-N 'Super- 




Figure 1. AzTEC and SCUBA coverage contours for the 
GOODS-N region demonstrates our uniform coverage. The dark 
rectangular contour corresponds to the AzTEC region with a map 
r.m.s. 5C1.16mJy at 1.1mm, the coverage region presented here. 
The grey contours, according to increasing line thickness, are the 
850 (im SCUBA contours for r.m.s. values of 4mjy, 2.5 mjy, and 
0.5 mjy respectively. The underlying map is the IRAC 3.6 fim im- 
age from the Spitzer legacy program (Dickinson et al. in prep). 
The AzTEC map represents a significant improvement in the uni- 
formity of coverage at faint flux levels. 



map.' The sensitivity variations of the AzTEC and SCUBA 
maps are compared in Fig. [T] 

In this work, we extract a catalogue of mm sources from 
the map and calculate number counts towards the faint end 
of the 1.1-mm galaxy population. The main results we dis- 
cuss here were obtained from the AzTEC data alone; data 
from other surveys have been used only as tools to check the 
quality of our map. A second paper will address counterpart 
identification of our AzTEC sources at other wavelengths 
(Chapin et al. in prep.). We present the JCMT/AzTEC ob- 
servations of GOODS-N in § [2] data reduction and analy- 
sis leading to source identification in § [3j properties of our 
source catalogue in § 31 the number counts analysis in § 
the discussion of results in § [6j and the conclusion in § [TJ 



2 AZTEC OBSERVATIONS OF GOODS-N 

AzTEC is a 144-element focal-plane bolometer array de- 
signed for use at the 50-m Large Millimetre Telescope 
(LMT) currently nearing completion on Cerro La Ne- 
gra, Mexico. Prior to permanent installation at the LMT, 
AzTEC was used on the JCMT between Nov. 2005 and 
Feb. 2006, p rimarily for deep , large-area blank field SMG 
surveys (e.g. IScott et alj|2008t Austermann et al. in prep.). 
We imaged the GOODS-N field at 1.1 mm with the AzTEC 
camera during this 2005-2006 JCMT observing campaign. 
Details of the AzTEC optical design, d etector array, and in - 
strument performance can be found in I Wilson et~al] (2008). 
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Each detector has a roughly Gaussian-shaped beam on the 
sky with an 18-arcsec full- width at half- maximum (FWHM). 
Given the beam separation of 22 arcsec, the hexagonal close- 
packed array subtends a "footprint" of 5 arcmin on the sky. 
Out of the full array complement of 144 bolometer-channels, 
107 were operational during this run. 

We mapped a 21 arcmin x 15 arcmin area centred on 
the GOODS-N field (12 h 37 m 00 s , +62°13'00") in unchopped 
raster-scan mode, where the primary mirror scans the sky 
at constant velocity, takes a small orthogonal step, then 
scans with the same speed in the opposite direction, re- 
peating until the entire area has been covered. We used a 
step size of 9 arcsec in order to uniformly Nyquist-sample 
the sky. We scanned at speeds in the range 60 arcsec s _1 - 
180 arcsec s _1 as allowed by the fast time constants of our 
micro-mesh bolometers, with no adverse vibrational system- 
atics. In total, we obtained 50 usable individual raster-scan 
observations, each taking 40 minutes (excluding calibration 
and pointing overheads). The zenith opacity at 225 GHz is 
monitored with the CSO tau meter, and ranged from 0.05- 
0.27 during the GOODS-N observations. This corresponds 
to 1.1 mm transmissions in the range 70-94%. A detailed de- 
scription an d justification of th e scan strategy we used can 
be found in I Wilson et""ai1 i|2008l '). 



3 DATA REDUCTION: FROM 

TIME-STREAMS TO SOURCE CATALOGUE 

In this section, we summarise the processing of the 
AzTEC/GOODS-N data, which is specifically geared to- 
wards finding mm point sources. The data red uction pro- 
cedure generally follows the method outlined in IScott et all 
(2008), although we emphasise several new pieces of analy- 
sis that were facilitated by the improved depth of this map 
over the COSMOS survey. We begin with the cleaning and 
calibration of the time-stream data in § 13.11 which includes a 
new investigation into the sample length over which to clean 
the data. In § 13.21 we describe the map making process and 
the optimal filtering for point sources. We asses the prop- 
erties and quality of the AzTEC/GOODS-N map in § 1331 
The depth of this survey has enabled us to ascertain the 
degree to which our data follow Gaussian statistics and de- 
tect, directly, a departure from it at long integration times 
indicating a component of signal variance due to source con- 
fusion. The astrometry of the map is analysed in § 13.41 and 
we describe the extraction of sources from the optimally fil- 
tered map in § 13.51 



3.1 



Filtering, cleaning, and calibration of 
time-stream data 



The AzTEC data for each raster-scan observation consists 
of pointing, housekeeping (internal thermometry, etc.), and 
bolometer time-stream signals. Because the bolometer data 
are sampled at 64 Hz, all other signals are interpolated to 
that frequency as needed by the analysis. The raw time- 
streams of the 107 working bolometers are fir st despiked 
and lo w-pass filtered at 16 Hz, as described in IScott et alj 
(2008). The despiked and filtered time-streams are next 
"cleaned" using a principal component analysis (PCA) ap- 
proach, which primarily removes the strong atmospheric sig- 



nal from the data. This "PC A-cleaning" method was devel- 
oped by the Bolocam group jLaurent et al. 20051 ) and later 
adapted for AzTEC, as described in lScott et all (|2008l ). As 
explained there, we also generate PCA-cleaned time streams 
corresponding to a simulated point source near the field cen- 
tre, in order to produce the point-source kernel, which is used 
later for beam-smoothing our maps (see § 13. 2f) . 

In this work we go beyond the analysis in IScott et al.l 
(2008) to verify that we have made good choices with re- 
gard to several aspects of the general cleaning procedure 
that has been adopted for all of the existing AzTEC data. 
We examine two outstanding questions in particular: 1) does 
PCA-cleaning work better than a simple common-mode sub- 
traction based only on the average signal measured by all 
detectors as a function of time? and 2) over what time scale 
should each eigenvector projection be calculated in order to 
give the best results? 

The first question addresses whether simple physical 
models may be used in place of PCA-cleaning, where the 
choice of which modes to remove from the data is not phys- 
ically motivated. We investigate this by creating a simple 
sky-signal template as the average of all of the detectors at 
each time sample. We then fit for an amplitude coefficient of 
the template to each detector by minimising the r.m.s. be- 
tween the scaled template and the actual data. This scaled 
template is removed from the bolometer data and we exam- 
ine the residual signal, which ideally consists only of astro- 
nomical signal and white noise. We find that this residual 
signal contains many smaller detector-detector correlations 
that are clearly visible in the data and are dominant com- 
pared to the signal produced by astronomical sources in the 
map. The residual time-stream r.m.s. from the simple sky- 
template subtraction is usually about twice the r.m.s. re- 
sulting from PCA cleaning. This test shows that the simple 
common- mode removal technique is insufficient. 

In the "stan dard" PCA-cleani ng procedure for AzTEC 
data, outlined in lScott et al.l (|200ST ) . the eigenvector decom- 
position is performed on each scan (~ 5-15 s of data). We 
now study which time scales give the best results using a 
statistical correlation analysis. We generate a bolometer- 
bolometer Pearson correlation matrix using sample lengths 
that range from a fraction of a second to tens of minutes 
(the length of a complete observation). On the shortest time 
scales, the correlation coefficients have large uncertainties 
due to sample variance (too few samples from which to make 
estimates). On time scales corresponding to a single raster- 
scan (~5-15 sec), however, the sample variance decreases 
and a clear pattern emerges: the strength of the correla- 
tions drops off uniformly with physical separation between 
the detectors. The most obvious trend is the gradient in 
correlations that we see with detector elevation, which is 
presumed to be produced by the underlying gradient in sky 
emission. As the sample length increases, a different pattern 
emerges, in which the dominant correlation appears to be 
related to the order in which the detectors are sampled by 
the read-out electronics, rather than their physical separa- 
tion. These correlations, likely due to electronics-related 1// 
drifts, are effectively removed when using scan-sized sample 
lengths (5-15 s) as well, since they appear as DC baseline 
differences on these short time-scales. These results verify 
that scan-sized sample lengths produce the best results as 
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they provide a sufficient number of samples on short enough 
time scales. 

After PCA-cleaning the bolometer signals, we apply a 
calibration factor to convert the bolometers' voltage time- 
streams int o units of Jy per bea m. Details of this procedure 
are given in I Wilson et alj (|200Sl ). The total error on the cal- 
ibrated signals (including the error on the absolute flux of 
Uranus) is 11%. 

3.2 Map-making and optimal Altering 

The map-making process used to generate the final opti- 
mally fi ltered AzTEC/GO ODS-N map is identical to that 
used in IScott et"al I (|2008l V and the reader is directed to 
that paper for the details of this process, which we briefly 
summarise below. 

We first generate maps for each of the 50 individual 
raster-scan observations separately by binning the time- 
stream data onto a 3" x 3" grid in RA-Dec which is tan- 
gent to the celestial sphere at (12 h 37 m 00 s , +62°13'00")- We 
chose the same tangent point and pixel size as t hat used for 
the S CUBA map of GOODS-N (see for example IPope et all 
l200rj ) so that the two maps can easily be compared in a fu- 
ture paper. We find that this pixel size provides a good com- 
promise between reducing computation time, while sampling 
with high resolution the 18-arcsec FWHM beams. Individual 
signal maps and their corresponding w eight maps for eac h 
observation are created as described in IScott et al.1 ll 20081) . 
along with kernel maps that reflect how a faint point source 
is affected by PCA-cleaning and other steps in the anal- 
ysis. Next, we form a single "co-added" signal map from 
the weighted average over all 50 individual observations. 
An averaged kernel map is also created in a similar way. 
The total weight map is calculated by summing the weights 
fro m individual obse rvations, pixel by pixel. As described 
bv lScott et al] (|2008l 1 we also generate 100 noise realization 
maps corresponding to the co-added map. 

We then use a spatial filter to beam-smooth our map 
using the point-source kernel, by optimally weighting each 
spatial-frequency component of this convolution accord- 
ing to the spatial power spectral density (PSD) of noise- 
realizati on maps. Details o f this optimal filter can also be 
found in lScott et all \200j \. 

3.3 Map quality: depth, uniformity, point-source 
response, and noise integration 

The final co-added, optimally filtered signal map for the 
GOODS-N field is shown in Fig.d Of the 315-arcmin 2 solid 
angle scanned by the telescope boresight during our survey, 
we expect ~250 arcmin 2 to be imaged uniformly by the com- 
plete AzTEC array. We identify this region by imposing a 
coverage cut. We find that weights within 70% of the cen- 
tral value occur in a contiguous region of 245 arcmin 2 . The 
map of Fig. [2] has been trimmed to only show this region. 
Much of the analysis presented here is limited to this re- 
gion, which we will henceforth refer to as the the "70% cov- 
erage region." The 1-a flux-density error estimates in the 
trimmed map range from 0.96 mjybeam -1 in the centre to 
1.16 mjybeam -1 at the edges. 

We also run the co- added kernel map through the same 
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Figure 3. Cross section of the point-source kernel. The Gaussian 
that best fits the inner R = 10 arcsec region is shown in the lighter 
shade and has a FWHM of 19.5 arcsec. The negative ring around 
the centre and other peripheral features (not visible here) are 
induced by PCA-cleaning as well as the optimal filter. 

filtering process as the signal map. The resulting filtered 
kernel map, whose profile is shown in Fig. [3] is our best 
approximation of the shape of a point source in the co- 
added, filtered signal map. As demonstrated in § 13.41 our 
pointing jitter/uncertainty has a sub-2-arcsec characteristic 
scale; this will have little impact on the kernel shape and 
therefore is not included in generating the kernel map. The 
negative troughs around the central peak are due to a combi- 
nation of array common-mode removal in the PCA-cleaning 
and de- weighting of longer spatial wavelength modes by the 
optimal filter. The point-source kernel also has radial scan- 
oriented features, or "spokes," due to PCA cleaning that 
are <0.1% of the kernel amplitude. The directions of these 
spokes would vary across the map as the scan angle changes 
with RA-Dec. Therefore, the kernel map accurately reflects 
these spokes only for point sources near the centre of the 
field. However, because it is difficult to analytically model 
a point source (through PCA cleaning and optimal filter- 
ing) and because the radial features are very faint, we use 
the kernel map as a point source-template for injection of 
sources in the simulations described later. 

Because this GOODS-N survey is the deepest blank- 
field survey conducted thus far with AzTEC on the JCMT, 
we demonstrate in Fig. [5] how the map noise averages down 
with the successive co-addition of individual observations. 
The central 200" x 200" region of the signal map and the 
noise realisation maps are used for this calculation. The x- 
axis represents the average weight of a 3-arcsec pixel in this 
region prior to filtering. A scale factor converts this raw 
weight to an effective time, T* , so that the final effective 
time equals the final integration time devoted to an average 
3-arcsec pixel in this central patch. Thus, the increment in 
T* gained with the addition of an individual observation is 
the effective integration time contributed by that particu- 
lar observation to the central region. The ith y-axis value is 
calculated by co-adding (averaging) individual signal maps 
from observations 1 through i, then applying the optimal 
filter, and finally taking the standard deviation of this co- 
added, filtered map in the central region. The crosses repre- 
sent the signal map. The 100 curves shown in a lighter shade 



An AzTEC 1.1 mm Survey of the GOODS-N Field I: Maps, Catalogue, and Source Statistics 5 



O 
O 

o 

r I 



EZ3 
CD 
CD 

tH 

CD 

O 
CD 



62.40 1 
62.35 I 
62.30 I 
62.25 ' 7 

L 

62.20 I 
62.15 I 
62.10 - 



© 9 O 23 



o 



020 

1* 



' J 34 
• /3i 



111 



1 O l6 © i5 



o 



8.0 



6.5 



5.0 



3.5 



2.0 



0.5 



-1.0 



189.4 189.2 189.0 

RA [degrees] J2000 



mJy 



Figure 2. AzTEC/GOODS-N signal map with the 36 S/NJJ3.5 source candidates circled. Information about these source candidates 
is given in Table [T] Here and in that table, source candidates are numbered in decreasing order of S/N. The source candidates marked 
with dashed-line circles do not belong to the robust sub-list, indicated by a horizontal line in Table [T] The map has been trimmed to 
show only the 70% coverage region (245 arcmin 2 ). 



are calculated by carrying out the same process on 100 noise 
realisations. In the absence of systematics or astronomical 
signal, we expect all curves to scale as 1/ \fT*, in accordance 
with Gaussian statistics, as indicated by the dashed line. At 
higher T* , we may expect a slight steepening in all curves 
because later co-additions better reflect our assumptions of 
circular symmetry (in the optimal filtering process) as we 
add more scan directions to the mix. However, this effect 
appears to be unmeasurably small in our data. 

While the noise realisations follow the 1/yT* trend, 
the signal map initially follows it but flattens near the point 
where ~20-30 s of effective time is spent on a 3-arcsec pixel. 
Switching the order in which signal maps are co-added does 
not alter this trend or the noisy behaviour of these points at 
large T* . Therefore, we conclude that: 1) single individual 
observations yield maps that are consistent with our noise 
realisations; 2) map features that do not survive scan-by- 
scan "jack-knifing," presumably astronomical signal due to 



source confusion, prevent the signal map's r.m.s. from im- 
proving as 1/VT*; and 3) the fact that noise realisations 
continue to follow this trend indicates that we are far from 
a systematics floor due to atmospheric or instrumental ef- 
fects, even at the highest T*. 



3.4 Astrometry calibration 

The pipeline used to produce this map of GOODS-N in- 
terpolates pointing offsets inferred from regular observa- 
tions of pointing calibrators interspersed with science tar- 
gets l|Wilson et al.ll2008l ; [Scott et aij|2008h . In order to ver- 
ify the quality of this pointing model for GOODS-N, both in 
an absolute sense, and in terms of small variations between 
passes, we compare the AzTEC m ap with the ex tremely 
deep 1.4 GHz VLA data in this field (|Richardsll200ol . Morri- 
son et al. in prep.). The radio da ta reduction and so urce list 
used here is the same as that of iPope et alj (|2006h . with a 
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Figure 4. Behaviour of the signal map's r.m.s. (crosses), as 
well as the r.m.s. of 100 separate noise realisations (collection 
of curves), as a function of the mean effective integration time 
T* spent on each 3-arcsec central pixel of map. The dashed curve 
shows the 1 / \JT* relationship expected in the absence of system- 
atics and astronomical signal. This demonstrates how the map 
noise averages down with the successive addition of more obser- 
vations. The "flattening" of the central r.m.s. at large T* in the 
signal map, compared to the noise maps, is due to astronomical 
signal. The fluctuations of this curve at large T* are simply due to 
noise in the r.m.s. itself, as re-ordering observations gives similar 
features near the same region. 



noise of ~5.3 /iJy at the phase centre. The catalogue is 
constructed with a 4-cr cut, and has positional uncertainties 
~ 0.2" (Morrison et al. in prep.). 

We stack the signal in the AzTEC map at the positions 
of radio sources to check for gross astrometric shifts in the 
AzTEC pointing model, as well as any broadening in the 
stacked signal which may indicate significant random offsets 
in the pointing between visits. A more detailed comparison 
between the mm and 1.4 GHz map is presented in (Chapin 
et al. in prep.) to assist with the MIR/NIR identifications of 
individual AzTEC SMGs, and the production of radio-NIR 
SEDs. 

The stack was made from the 453 1.4 GHz source posi- 
tions that ar e within the unifor m noise region of the AzTEC 
map. As in IScott et all (|2008h we check for an astromet- 
ric shift and broadening by fitting a simple model to the 
stacked image, which consists of an astrometric shift (<5RA, 
<5Dec) to the ideal point source kernel, convolved with a 
symmetric Gaussian with standard deviation a p . This Gaus- 
sian represents our model for the random pointing error in 
the AzTEC map. We determine maximum likelihood esti- 
mates SRA = 0.2", <5Dec = -0.9", and cr p = 0.6". The 
expected positional uncertainty (in each coordinate) for a 
point source with a purely Gaussian bea m is approximately 
0.6xFWHM/(S/N) (see the Appendix in llvison et aD^OOTt ) 
where the FWHM is 18arcsec in our case. The S/N of our 
stack is approximately 10, so the expected positional uncer- 
tainty is ~ 1". Therefore the total astrometric shift measured 
by the fitting process, 0.9", is consistent with the hypothe- 
sis that there is no significant underlying shift. We also note 
that the \ 2 function for this fit is extremely shallow along 
the cr p axis, so although the minimum occurs at 0.6", it is 
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Figure 5. Pixel flux histogram of the final signal map in a dark 
shade and the average pixel flux histogram made from 100 noise 
realizations in a lighter shade. The positive tail and smaller neg- 
ative excess in the signal map is due to the presence of point 
sources. 



not significantly more likely than 0". We therefore conclude 
from this analysis that there is no significant offset, nor beam 
broadening caused by errors in the pointing model. 



3.5 Source finding 

To investigate the presence of astronomical sources in our 
map, we plot in Fig.[5]a histogram of pixel fluxes in the 70% 
coverage region of the field. Also shown in a lighter shade 
is the average pixel histogram made from the 100 noise- 
realization maps. The noise histogram can be modelled well 
by a Gaussian centred on OmJy with a standard deviation 
of 1.0 mjy. The obvious excess of large positive pixel values 
and the small excess of negative values in the signal map are 
caused by the presence of sources. 

To identify individual point sources, we first form a S/N 
map by multiplying the final (i.e. co-added and filtered) sig- 
nal map by the square-root of the weight map. We then iden- 
tify local maxima in this S/N map with S/N 3.5. There 
are 36 local maxima that meet this condition in the 70% 
coverage region of the field. Our analysis of these source 
candidates is simplified because no pair of them are close 
enough to significantly alter each other's recovered flux den- 
sities (>36 arcsec apart in each case). We have evidence that 
AzGNOl is a blend of two sources. However, since this knowl- 
edge is not based on AzTEC data alone, we defer a detailed 
discussion of that source for the second paper of this series 
(Chapin et al. in prep.). 

The final signal map and these source candidates 
are shown in Fig. Table Q] lists details of all the 
AzTEC/GOODS-N ^ 3.5-ct source candidates, including 
their locations, measured fluxes, S/N, and additional quan- 
tities which are defined below. The source positions are 
given to sub-pixel resolution by calculating a centroid for 
each local maximum based on nearby pixel fluxes. Sources 
with clear counterparts in the SCU BA map of GOODS-N 
l|Borvs et al.ll2003l : IPope et alJl2005l ) are highlighted in Ta- 
ble en 
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Table 1. Source candidates in AzTEC/GOODS-N with S/N>3.5 ordered according to S/N. The horizontal line between AzGN28 and 
AzGN29 represents our threshold for source robustness, as explained in The last two columns are defined in § 14,21 The superscripts 

S and N highlight sources in our robust sub-list that lie within the considered SCUBA region (where the 850-/jm r.m.s. is <2.5mjy). 
The sources denoted by S have robust detections at 850 (im within 12 arcsec of the given positions while the sources denoted by N do 
not (Chapin et al. in prep.). 



4 THE AZTEC/GOODS-N SOURCE 
CATALOGUE 

As evident from Table [1] the number of source candidates 
increases rapidly with decreasing S/N. However, if we use 
a S/N threshold to make a sub- list of the sources in Ta- 
ble [1] the false positives contained in such a list will also 
increase with lower S/N thresholds. Our aim here is to find 
a S/N threshold above which >95% of source candidates are, 
on average, expected to be true sources. This is a practical 
choice aimed at maximising the number of sources recom- 
mended for follow-up studies (the subject of Chapin et al. 
in prep.) in a way that limits the effect of false detections 
on any conclusions drawn. The horizontal line in Table [1] 
below source AzGN28 (S/N^3.75) marks the cut-off of the 
sub- list that we expect will satisfy our robustness condition. 



We first explain in § 14. li the analysis of false detection rates 
(FDRs) that yields this threshold. In that sec tion, we go 
beyo nd previous FDR treatments for AzTEC (|Scott et al.l 
2008) and derive some general results about FDRs that are 
applicable to (sub)mm surveys in general. 



Next, we explain in § !4.2l the last two columns of Table[T] 
which contain a re-evaluation of source flux densities and an 
assessment of the relative robustness of our source candi- 
dates. Then, in § 14.31 we discuss the survey completeness 
and present a brief consistency check of our source candi- 
dates against SCUBA detections at 850 /J,m. 
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Source Threshold 
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Pure-noise Negative FDR 


4.55 


1.58 


0.33 


0.00 


Best-fit-model Negative FDR 


5.96 


2.85 


1.16 


0.04 



Table 2. The number of source candidates passing a given S/N 
threshold in the actual map are indicated in row 1. Several meth- 
ods for determining the false detection rates (FDRs) were ex- 
plored. "Pure-noise" refers to averages computed over 100 noise- 
realization maps. "Best-fit-model" corresponds to averages from 
100 noisc+source realization maps using the best fit model of 
§ 15,11 We have settled on the values of row 2 as our nominal FDRs 
because they give a conservative overestimate, as explained in the 
text. 



4.1 False detection rates 

Two obvious methods for estimating the false detection rate 
(FDR) of a survey are to run the source finding algorithm 
on: 1) simulated noise realization maps; or 2) the negative of 
the observed signal map. For several S/N thresholds, Tabled 
lists the number of source candidates in the actual map (row 
1), the average number of "sources" found in simulated pure 
noise realizations (row 2), and the number of "sources" in 
the negative of the actual map (row 4). When using the map 
negative, regions within 36 arcsec of a bright positive source 
were excluded in order to avoid their "negative ring" (see 
Fig. El). 

We conclude that these two estimates of the FDR are 
not very accurate for our maps. Because of the high number 
density of SMGs in the sky compared to our beam size, ev- 
ery point of the map is in general affected by the presence of 
sources. This source confusion causes the simple FDR esti- 
mates above to be inaccurate. In particular, there are equal 
numbers of negative and positive "detections" in noise real- 
izations to within the statistical error of our noise simula- 
tions, as indicated by rows 2 and 5. However, the presence 
of sources skews this balance in the actual map, making the 
false negatives rate higher than the pure-noise numbers and 
the false positives rate (what we are after) lower than the 
pure-noise numbers. 

Both these effects can be understood by considering 
the following hypothetical construction: a noise-less AzTEC 
map of the sky containing many point sources, all with the 
shape of the point-source kernel. Because each kernel has 
a mean of zero, such a map would have an excess of neg- 
ative valued pixels over positive valued pixels (about 70% 
to 30%) to counter the high positive values near the centre 
of the kernel (see [3}. When noise that is symmetric around 
zero is "added" to such a map, this small negative bias will 
cause a larger number of high-significance negative excur- 
sions in that sky map compared to a map containing just 
the symmetric noise. The pixel flux histogram of the actual 
map, shown in Fig. [5] (darker shade), also shows evidence 
of this effect through its negatively shifted peak as well as 
the excess of negative pixels in comparison with pure noise 
realizations (lighter-shade histogram). This small negative 
bias, in pixels that do not lie atop a source peak, also ex- 



plains why there are fewer high-significance false positives 
in an actual sky map compared to a pure noise map. 

To verify our reasoning, we generated 100 noise+source 
realizations for the best-fit number counts model described 
in § 15.11 For each realization, we find the number of posi- 
tive and negative "detections" just as for the true map. False 
positives are defined as detections occurring >10 arcsec away 
from inserted sources of brightness >0.1mJy. The FDR re- 
sults for these simulations are given in rows 3 and 6 of Ta- 
ble [2] The results show that the negatives rate is indeed 
boosted by the presence of sources, compared to pure noise 
maps (rows 2 and 5). Furthermore, the negative FDR of the 
actual map (row 4), which drops to at a S/N of 4.2, is sta- 
tistically consistent with the simulated negative FDR means 
of row 6. As expected, the simulated false positives rate is 
lower than the pure-noise FDR, as evident from row 3. 

As the true positive FDR depends on the number 
counts, we adopt the model-independent pure-noise values of 
row 2 as our nominal FDRs. These will be conservative over- 
estimates of the FDR regardless of the true 1.1 mm number- 
counts of the GOODS-N field. 

Based on these nominal FDRs, we divide the source can- 
didate list of table [1] into two categories of robustness, with 
the dividing line at a S/N of 3.75. On average, we expect 
1-2 source candidates with S/N^3.75 (above the horizontal 
line in Table [TJ and 1-3 candidates with S/N<3.75 (below 
the line) to be false detections. 

4.2 Flux bias correction 

In our map, where the signal from sources does not com- 
pletely dominate over noise, the measured flux density can 
be significantly shifted from the true 1.1mm flux density 
of a source due to noise. The measured flux densities in 
column 4 of Table Q] are more likely to be overestimates 
than underestimates of the true flux densities because of the 
sharply decreasing surface density of (sub)mm galaxies with 
increasing flux density. As th is slope in the number counts is 
quite steep (see for example Blain ct al. 1999; Bargcr "et al.l 
19991: lEales et alJl200Cl : iBorvs et al.ll2003l : iGreve et alj|2004l ; 
Coppin et al.ll2006h . this bias can be a large effect. There- 
fore, we estimate a "de-boosted" flux density for all our 3.5-<r 
source candidates. This estimate i s base d on the Bayesian 
technique laid out in ICoppin et all (2005) for calculating the 
posterior flux density (PFD) distribution of each source. 

The number-counts model that we use to generate the 
prior is given by 



(1) 



dS S 

where dN/dS represents the differential number counts of 
sources with flux density 5*. We use N' = SSOOmJy-Meg -1 
and S' — 1.5 mjy, which is consist ent with taking the 
Schechter function number-counts fit of lCoppin" et al.l (2006) 
and scaling the 850 /im fluxes by a factor of 2.2 to approx- 
imate the 1.1mm fluxes of the same population. It is suffi- 
cient to use a prior that is only approximately correct, since 
many of the derived results (as we have checked explicitly) 
are independent of the exact form of the assumed number 
counts. We take as our Bayesian prior the noise-less pixel 
flux histogram of a large patch of sky simulated according 
to this model. Since our point-source kernel has a mean of 
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Figure 6. Survey completeness for the S/NJJ3.75 cut used here 
to select robust sources is represented with the dark symbols and 
error bars. The lighter symbols and error bars are estimates of 
the survey completeness when the integrated posterior flux dis- 
tribution below Omjy is required to be <5%. 



zero, the prior is non-zero for negative fluxes and peaks near 
OmJy. 

The de-boosted flux density given in column 6 of Table[T] 
is the location of the PFD's local maximum closest to the 
measured flux density. This de-boosted flux density is fairly 
insensitive to changes in the prior that correspond to other 
number-counts models allowed by current constraints. The 
upper and lower error bounds quoted for a de-boosted flux 
density correspond to the narrowest PFD interval bracketing 
the local maximum that integrates to 68.3%. 

In order to determine the relative robustness of each 
source individually, we calculate the integral of the PFD be- 
low zero flux. This quantity, given in column 7 of Table [T] 
is not a function of just S /N but depends on the flux (sig- 
nal) and its error (noise) separately. Although the values 
given in column 7 can vary appreciably among reasonable 
choices of number-counts priors and the PFD integration 
upper bounds (set to zero here), the source robustness order 
inferred by the non-positive PFD integral is quite insensi- 
tive to these choices. Therefore, the values in column 7 pro- 
vide a useful indicator of the relative reliability of individual 
sources. 

However, due to the arbitrariness present, the values in 
this column cannot be used to directly calculate the FDR of 
a source list. For instance, the sum of column 7 values for 
our robust source list is ~0.5, which is an underestimate of 
the expected FDR (see 8 14.10 . We note that, for our choice of 
prior, the requirement of a non-positive PFD ^5% happens 
to identify the same robust source-candidate list as the S/N 
cut of 3.75. However, this statement is specific to a particular 
choice of prior and PFD integration upper bound. 



4.3 Survey completeness and comparison with 
SCUBA detections 

We next compute the survey completeness by injecting one 
source at a time, in the form of the point-source kernel scaled 
to represent each flux, at random positions in the GOODS- 
N signal map (Fig. [2} and tallying the instances when a new 



source is recovered with S/NJ53.75 within lOarcsec of the 
insertion point. We choose this radius because it is small 
enough for conducting quick searches in our simulations 
and because, barring incompleteness, simulations show that 
>99.5% of ^ 3.75-rj sources will be found within lOarcsec 
of their true position given the size of the AzTEC beam and 
the depth of coverage. This method of calculating complete- 
ness allows for the inclusion of "confusion noise" without 
altering the map properties appreciably, because only one 
artificial source i s injected per simulation l|Scott et alj 20061 ; 
IScott et al]|2008h . 

We have also assessed completeness by inserting point 
sources of known flux, one at a time, into pure noise- 
realisation maps rather than the signal map. With this 
method, we also require that each recovered artificial source 
has a < 5% non-positive PFD. Since this constraint is es- 
sentially equivalent to a limiting S/N threshold of 3.75 (as 
evident from Table [TJ , it is not surprising that the survey 
completeness determined this way (lighter-shade points of 
Fig. [6} is similar to that derived from the previous method. 
The similarity in results also shows that the effects of con- 
fusion noise on survey completeness is small. 

Finally, to verify that our source-candidate list has over- 
lap with previously detected extragalactic (sub)mm sources, 
we compare our source list against 850-/im SCUBA detec- 
tions within overlapping survey regions. For this purpose, 
we only consider the regions in the SCUBA/850-/im map 
with noise r.m.s. < 2.5 mjy. Of the 28 AzTEC sources in 
the robust list, 11 lie within this region of the SCUBA 
map; of these 8 (73%) have robust detections at 850 /im 
l|Pope et al.ll2005h within 12arcsec of the AzTEC position. 
Those 8 are highlighted with the superscript "S" in Table [1] 
while the other 3 are marked with the superscript "AT." On 
the other hand, all 38 robust SCUBA sources within the 
r.m.s. < 2.5 mjy region |Pope et al.ll2005l ; IWall et al.ll2008h 
lie within the 70% coverage region of AzTEC. In Chapin 
et al. (in prep.) we will discuss the 850 /im properties of 
AzTEC sources by performing photometry in the SCUBA 
map at AzTEC positions, and more fully explore the overlap 
of the AzTEC and SCUBA populations in general. 



5 1.1 MM NUMBER COUNTS 

Using our AzTEC/GOODS-N data, we next quantify the 
number density of sources as a function of their intrinsic 
(de-boosted) 1.1mm flux. These counts cannot be read di- 
rectly from the recovered distribution of source flux densities 
due to: 1) the bias towards higher fluxes in the data (as de- 
scribed in section [472]) . which includes false detections; and 
2) the survey incompleteness at lower fluxes. In order to esti- 
mate the counts we use two independent methods: a Monte 
Carlo technique that implicitly includes the flux bias and 
completeness issues; and a Bayesian approach that accounts 
for both these effects explicitly. 

Fig. [7] shows the results of our number-counts simula- 
tions. It shows the source flux density histogram simulated 
for the best fit model from the parametric method overlaid 
on the actual distribution from the true map. It also shows 
the differential number counts vs. de-boosted source flux 
density as returned by both methods. The dot-dashed lines 
in the lower right correspond to the survey limits of the 
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frequentist and Bayesian approaches, which are 27.8 and 
33.8 deg -2 mjy -1 , respectively. The survey limit is the y- 
axis value (number counts) that experiences Poisson devia- 
tions to zero sources per mjy-bin 32.7% of the time, given 
the map area considered. The two limits differ slightly be- 
cause the frequentist simulations include the slightly larger 
area 50% coverage region, as opposed to the 70% coverage 
region that we use for the Bayesian method. The survey 
limit occurs at around 6 mjy for both the best-fit frequentist 
and Bayesian type simulations. Thus, we are not sensitive to 
the differential number counts with 1 mjy resolution beyond 
that point. 

The power of the AzTEC/GOODS-N survey is in con- 
straining number-counts at lower flux densities, given the 
depth reached in this relatively small field. We have, how- 
ever, excluded results below the <2mJy level from both 
methods, because of low survey completeness (<10%) and 
the possibility of increasing systematic effects. Therefore, 
the noteworthy features of Fig. are the points from the 
Bayesian approach, indicated by crosses and error bars, in 
the range 2 mjy to 6 mjy and the allowed functional forms 
from the parametric (frequentist) method within those flux 
density bounds. Models allowed by the 68.3% confidence 
interval of the parametric method form the shaded region 
while the dark curve is the best-fit model. Given the error 
bounds from the two methods, they are in good agreement. 
Both methods are briefly described below. 

5.1 Parametric frequentist approach 

An obvious choice of indicator for the underlying source pop- 
ulation is the recovered brightness distribution of source 
candidates in the GOODS-N map. Here, we use a S/N 
threshold of 3.5 and the 50% coverage region of the map. 
After identifying S/N^3.5 source candidates, we make a 
histogram of their measured flux densities using 0.25 mjy 
bins, for comparison against histograms made from simulat- 
ing various number-counts models. Th is approach is s i milar , 
in spirit, to the method employed in lLaurent et al.l (|2005l ) 
and the parametric version of number counts derived in 
ICoppin et af] |2006). However, we avoid intermediate ana- 
lytical constructs, as the procedure outlined below accounts 
for all relevant effects. 

We generate model realisation maps by injecting kernel- 
shaped point sources into noise realisation maps. The input 
source positions are uniformly distributed over the noise re- 
alisation map while their number density and flux distri- 
bution reflect the number-counts model being considered. 
For every model we have considered, we make 1200 simu- 
lated maps by constructing 12 different source realisations 
for each of the 100 noise realisation maps. Next, we use the 
same source-finding algorithm used on the signal map to ex- 
tract all S/N^3.5 peaks in each simulated map. We then 
compare the average histogram of recovered source fluxes 
from the 1200 model realisations against the actual distri- 
bution of source fluxes. The data vs. models comparison 
is restricted to the 3.5-8 mjy measured flux density range. 
This comparison process is illustrated in Fig. [7] 

The likelihood of the data given a model is dete r mined 
acco rding to Poisson sta tistics as in lLaurent et al.l (|2005l ) 
and ICoppin et al.l |2006). One set of parameterised mod- 
els that we have explored has the functional form given by 




Source Flux at 1 . 1 mm [mjy] 

Figure 7. The thick solid curve and the enveloping shaded region 
correspond to the best fit number counts model and the 68.3% 
confidence interval from the parametric approach of § 15.11 The 
distribution of measured fluxes of 3.5-a sources in the actual map 
is shown by the triangles in the 3.5-8 mjy interval while the cor- 
responding average distribution of the best fit model is indicated 
by the thin solid-line histogram. The difference between the thick 
solid line and the thin solid histogram indicates the importance of 
accounting for flux boosting and completeness. The crosses and 
error bars represent the differential number counts derived from 
the Bayesian method, which are in excellent agreement with the 
result from the parametric method. The dashed-line curve indi- 
cates the Bayesian prior. The upper and lower dot-dashed lines 
indicate the survey limits of the Bayesian and parametric meth- 
ods, respectively. 



Equation[T] We chose to re-parametrise these models so that 
the normalisation factor depends on only one of the fit pa- 
rameters. The parameters we chose are the same S' as in 
EquationQ]and A^3 m j y , the differential counts at 3 mjy, given 
by 



N: 



3m Jy 



N' 



S' 
3mJy 



-3mjy / S 



In terms of these parameters, Equation [T] becomes 

*>L-N* T ( 3mJ y \ ,-(g-3mJy)/g' 

ds ~ 3mJy V S J 



(2) 



(3) 



We explored the S'-Ns^jy parameter space over 
the rectangular region bracketed by 0.5-2 mjy and 60- 
960 mjy" 1 degree -2 using a (AS', AA^mjy) cell size of 
(0.15,60). The likelihood function, C, is a maximum for 
the model with S' = 1.25 ± 0.38 mjy and Ar 3m j y = 300 ± 
90 mjy -1 de gree . We did not assume % -like behaviour of 
— ln(£) for calculating the 68.3% confidence contours whose 
projecti ons are the e r ror ba rs quoted above. Instead, as out- 
lined in iPress et al.l l|l992f ). we made many realisations of 
the best-fit model and put them through the same param- 
eter estimation procedure that was applied to the actual 
data. In terms of the goodness of fit, we find that 66% of 
the simulated fits yield a higher value of — ln(£) compared 
to the actual value. Fig.[7]shows this best-fit number-counts 
estimate against the de-boosted 1.1mm flux density along 
with a continuum of curves allowed by the 68.3% confidence 
region. 



An AzTEC 1.1mm Survey of the GOODS-N Field I: Maps, Catalogue, and Source Statistics 11 



(2004) 
et al. (2005) 
et ol. (2006] scaled 

rk (parameiric) 

5d fit 




1.1 mm Flux Density [mlyj 

Figure 8. The cumulative (integral) number counts from 
other 1.1-1.2 mm surveys are shown alongside our results. The 
AzTEC/GOODS-N parametric number-counts results are indi- 
cated by the hatched region that represents the 68.3% confi- 
dence region for parametric models. The dot-dashed line indi- 
cates the survey limit. Results from the Bolocam 1.1mm Lock- 
man Hole survey are indicated by a thin solid line and two 
bounding dotted-lines that repr esent the best-fit mode l and 68.3% 
confidence region as found by IMalonev etaL I j2005l) The 1.2- 
mm MAMBO-IRAM results reported in lGreve etTiTl|2004h also 
sh own (triangles ) . The stars represent the "reduction D" results 
of ICoppin et al with 850 fitn flux densities scaled by the 

factor 1/2.08 as explained in the text. The dashed curve indicates 
the best combined fit to the Bayesian results from both surveys. 



5.2 Bayesian method 

We also estimate number counts from the individual source 
PFDs calculated in § 14.21 using a modified vers i on of the 
bootstrapping method described in ICoppin et al] (|2006h . A 
complete discussion of the modifications and tests of the 
method will be presented in Austermann et al. (in prep.). 
For these calculations, we use only the sub-list of robust 
sources in Table Q] We have repeated this bootstrapping 
process 20,000 times to measure the mean and uncertainty 
distributions of source counts in this field. The differential 
and integrated number counts extracted with this method, 
using 1 mjy bins, are shown in Fig. [7] Our simulations show 
that the extracted number counts are quite reliable for a 
wide range of source populations and only weakly dependent 
on the assumed population used to generate the Bayesian 
prior (the dashed-line curve of Fig. with the exception of 
the lowest flux density bins, below 2 mjy, which suffer from 
source confusion and low (and poorly constrained) complete- 
ness. Overall, the results from the Bayesian method are in 
excellent agreement with those from the parametric method 
between the lower sensitivity bound (2mJy) and the survey 
limit (~6mJy). 



6 DISCUSSION 

In Fig. [8] we display our cumulative number-counts results 
with the 68.3%-allowed hatched region derived from the 



Survey 


S' 


^3mJy 


"dust 


AzTEC/GOODS-N 


1.25 ± 0.38 


300 ± 90 




AzTEC/GOODS-N 








+ SCUBA/SHADES 


1.60 ±0.25 


274 ± 54 


2.84 ± 0.32 




8' 


^3mJy 


"dust 


S' 


1 


0.05 


-0.32 


^3mJy 


0.05 


1 


-0.8 


a dust 


-0.32 


-0.8 


1 



Table 3. Best-fit Schechter function parameters and dust 
emissivity spectral index using the Bayesian results from the 
AzTEC/GOODS-N, SCUBA/SHADES, and combined surveys. 
The correlation matrix for the combined fit is also listed. Caveats 
on this analysis are given in the text. 



parametric method. We next compare those results with pre- 
vious surveys of (sub)mm galaxies. Combined results from 
the 1.2-mm MA MBO surveys of t he Lockman Hole and 
ELAIS-N2 region (|Greve et al]|2004l) and the 1.1 -mm Bolo- 
cam Lockman Hole survey (IMalonev et al]|2005t ) are shown 
in Fig. [8] Our GOODS-N number counts are in good agree- 
ment with MAMB O results. Our results are in disagreement 
with the results of lMalonev eTail (|2005l ). even within a lim- 
ited flux range such as 3-6 mjy where we expect both sur- 
veys to be sensitive to the number counts. 

In Fig. [5] we also compare our resu lts with the 850-/im 
number counts of ICoppin et al] ((2006) . If the 1.1-1.2 mm 
surveys detect the same population of sub-mm sources seen 
by SCUBA at 850 /im - an assumption that is not obviously 
valid given the po ssible redshift-dependent selection effects 
l|Blain et alj|2002j ) - we would expect a general correspon- 
dence between number counts at these two wavelengths, 
with a scaling in flux density that represents the spectral 
factor for an average source. Therefore, we perform a simul- 
taneous fit to the SCUBA/SHADES and AzTEC/GOODS- 
N differential Bayesian number counts in order to deter- 
mine the average dust emissivity spectral index, adust (and 
thus the flux density scaling factor from 1.1 mm wavelength 
to 850pm wavelength), and the parameters, A r 3 m j y and S', 
of Equation [3] This fit results in the best-fit parameters 
and correlation ma trix given in Table [3] We overlay the 
ICoppin et ail ([2006) number counts on Fig. [8] with the 850- 
/im fluxes scaled by the scaling factor derived from this fit, 
which is 2.08±0.18. For visual comparison, the shaded re- 
gion of Fig. |8j which represents our parametric result, is 
sufficient because it represents well the results from both 
methods (see Fig. 0. Fig. |S] shows that the scaled SCUBA- 
SHADES points fall well within the bounds allowed by our 
results. 

The adust of Table [3] was computed for the nominal 
AzTEC and SCUBA band centres, which are 1.1mm and 
850 /im respectively. However, the quoted error on «dust 
brackets the effects of small shifts in the effective band cen- 
tres due to spectral index differences between SMGs and flux 
calibrators. The dust emissivity spectral index may also be 
estimated by averaging the 1.1mm to 850 ^m flux density 
ratio of individual sources or by performing the appropriate 
stacking analysis. Due to the moderate S/N of sources in 
our surveys, the effects of flux bias and survey completeness 
must be accounted for in such analyses. Therefore, perform- 
ing a combined fit to the differential number counts vs. de- 
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boosted flux from the two surveys, where those effects are 
already included, is an appropriate method for estimating 
the spectral index. From Fig. [8] the hypothesis that SCUBA 
and AzTEC detect the same underlying source population 
appears plausible. 

However, we do not comment on the formal goodness 
of fit as the \ 2 obtained for the combined fit is unreason- 
ably small because the full degree of correlation between 
data points is underestimated i n the standard comp utation 
of the two covariance matrices (jCoppin et al.l [20061 ). In ad- 
dition, the best-fit parameters of the combined fit may have 
a large scatter from a global mean value (if one exists), due 
to sample variance, as the two surveys cover different fields. 
Although S CUBA 850 /xm nu mber-counts are available for 
GOODS-N (jBorvs et alJl2003h . the survey region (see Fig.[jjj 
and the method used to estimate number counts in that work 
are quite different from those used here. There fore, we chose 
to fit to the SCUBA/SHADES number counts (jCoppin et all 
2006) instead, since they were determined using methods 
similar to ours. 



ber counts alone cannot really test this hypothesis. A more 
thorough study of whether AzTEC is selecting a systemati- 
cally different population than SCUBA can come only from 
comparison of the redshifts and multi-wavelength SEDs of 
the identified galaxies, which we will describe in Chapin et 
al. (in prep.), the second paper in this series. 

There is also a survey of GOODS-N with MAMBO at 
1.25mm performed by Greve et al. (in prep.). A compari- 
son between these two millimetre maps and, possibly, the 
SCUBA 'Super-map' is reserved for a future paper (Pope et 
al. in prep.). 

This AzTEC/GOODS-N map is one of the large blank- 
field SMG surveys at 1.1 mm taken at the J CMT. Com- 
bined with the AzTEC surveys in the COSMOS (|Scott et all 
2008) and SHADES (Austermann et al. in prep.) fields, these 
GOODS-N data will allow a study of clustering and cosmic 
variance on larger spatial scales than any existing (sub)mm 
extragalactic surveys. 
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7 CONCLUSION 

We have used the AzTEC instrument on the JCMT to image 
the GOODS-N field at 1.1 mm. The map has nearly uniform 
noise of 0.96-1.16 mjy beam -1 across a field of 245 arcmin 2 . 
A stacking analysis of the map flux at known radio source 
locations shows that any systematic pointing error for the 
map is smaller than 1 arcsec in both RA and Dec. Thus, 
the dominant astrometric errors for the 36 source-candidates 
with S/N^3.5 are due to noise in the centroid determination 
for each source. Using a S/N^3.75 threshold for source ro- 
bustness, we identify a subset of 28 source candidates among 
which we only expect 1-2 noise-induced spurious detections. 
Furthermore, of the 11 AzTEC sources that fall within the 
considered region of the SCUBA/850-^im, 8 are detected un- 
ambiguously. 

This AzTEC map of GOODS-N represents one of the 
largest, deepest mm-wavelength surveys taken to date and 
provides new constraints on the number counts at the faint 
end (down to ~2mJy) of the 1.1 mm galaxy population. 
We compare two very different techniques to estimate the 
number density of sources as a function of their intrinsic 
flux-a frequentist technique based on the flux histogram 
of detected sources in the map similar in spirit to that of 
lLaurent et alj J2005I). and a Bayesian approach similar to 
that of ICoppin et all l|2006h . Reassuringly, the two tech- 
niques give similar estimates for the number counts. Those 
results are in good agr e ement with the number counts es- 

but differ significantly from 



timates of I Greve et al 



those of iMalonev et al 



The 1.1 mm number counts from this field are consistent 
with a direct fl ux scaling of the 850 [im SCUBA/SHADES 
number counts l|Coppin et al.ll2006l ) within the uncertainty 
of the two measurements, with a flux density scaling fac- 
tor of 2.08 ± 0.18. If we assume that the two instruments 
are detecting the same population of sources, we obtain a 
grey body emissivity index of 2.84 ±0.32 for the dust in the 
sources. While there is no evidence based on the number 
counts that 1.1 mm surveys select a significantly different 
population than 850 u.m surveys, we caution that the num- 
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